Confined active nematic flow in cylindrical capillaries 
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We use numerical modelling to study the flow patterns of an active nematic confined in a cylin- 
drical capillary, considering both planar and homeotropic boundary conditions. We find that active 
flow emerges not only along the capillary axis but also within the plane of the capillary, where 
radial vortices are formed. If topological defects are imposed by the boundary conditions, they act 
as local pumps driving the flow. At higher activity we demonstate escape of the active defects and 
flow into the third dimension, indicating the importance of dimensionality in active materials. We 
argue that measuring the magnitude of the active flow as a function of the capillary radius allows 
determination of a value for the activity coefficient. 

PACS numbers: 83.50.-v,47.50.-d,87.10.-e 



There is currently widespread interest in understand- 
ing active systems which are key to the physics of life 
processes ^\^. Active systems are materials that pro- 
duce their own energy, or respond to a continuous input 
of energy, so that they naturally operate out of thermo- 
dynamic equilibrium [5|. Examples range from bacte- 
rial suspensions [6(1, solutions of bio-filaments driven by 
molecular motors [3, active colloids [sl, [9[ , and cellular 
pattern formation jioj , to floc king birds and fish [ll|, [l2j 
and vibrated granular matter [13|. Despite the diversity 
of these systems there is recurring behaviour, such as pat- 
tern formation 13, [HI , the prevalence of instabilities [HI 



An im- 



and the appearance of topological defects [17|, 
portant question in formulating theories of active matter 
is to identify the extent to which these behaviours are 
generic and to then identify the dominant control pa- 
rameters. 

Approaches that have proved very useful in categoris- 
ing passive soft matter systems are continuum theories 
which rely primarily on the symmetries of the prob- 
al continuum theories have been pro- 



lems 

posed for active systems but matching theories to exper- 
iment, and giving a microscopic interpretation of the pa- 
rameters, is still in its infancy [2]. Noting that nematic 
ordering and topological defects have been observed in 
active systems [sl, [23|, [l^, one of the most interesting 
candidates for a continuum approach is closely related 
to the Beris-Edwards equations of liquid crystal hydro- 
dynamics H, Hll. The important addition, representing 
the activity, is a term in the stress proportional to the 
nematic tensor order parameter. Hence any gradient in 
the order parameter field initiates a flow. 

In the bulk, the activity drives instabilities which de- 
stroy the nematic order 19|, 25-23- However, confining 
surfaces can stabilise a nematic state and correspond- 
ing time-independent flow fields. This has been demon- 
strated in a flat cell geometry 28|, , in cells with static 
and moving boundaries ^o|. Other works have consid- 
ered defect driven hydrodynamics in active circular re- 
gions surrounded by a passive background and spon- 



taneous circulation in cytoplasmic streming [32|. 

In this Letter we consider confined active nematic flows 
in a cylindrical capillary, as a common example of an 
experimetnally relevant geometry. This is a very rich 
system as topological defects can emerge naturally, as a 
consequence of the boundary conditions. Solving the con- 
tinuum equations of motion in the flow aligning extensile 
regime we find that active flow emerges not only along 
the capillary axis but also within the plane of the capil- 
lary, where radial vortices are formed. Interplay between 
the primary and secondary velocities leads to a spiralling 
flow field with a well-defined pitch moving down the cap- 
illary. We show that topological defects, imposed by het- 
erogeneous boundary conditions, act as sources driving 
the flow. We also argue that measuring the magnitude 
of the active flow as a function of the capillary radius al- 
lows determination of a value for the activity coefficient. 

Active nematic flow is approached by using numerical 
modelling based on continuum equations for apolar active 
materials 
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29|. The model characterises the orienta- 
tional order of active fluids by a symmetric and traceless 
order parameter tensor Q^j, with the largest eigenpair 
representing the magnitude and direction of the local ori- 
entational order. In parallel, the flow of the active fluid 
is described by a velocity flow field ui and density p. The 
dynamic equations for Qij and Ui constitute the two gov- 
erning equations of motion for the active nematic: 



{dt + Ukdk)Qij 

p{dt + Ukdk)ui 



Sij — 



VHij 



(1 



rjdj [diUj 
3dpP)dkUkSij], 



(1) 

(2) 



where dt is a derivative in time t, 9^ is a derivative in the 
Cartesian spatial coordinate x, ^, F is a rotational dif- 
fusion coefficient, r] is the viscosity, and A is a first active 
parameter which can be considered just as a renormali- 
sation of the free energy. Summation over repeated in- 
dices is assumed. The leading contribution to the active 
dynamics is controlled by the parameter Q in the active 
contribution to the non-Newtonian terms of the stress 
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tensor 
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QikHkj — HikQkj^ (3) 
(4) 



where ^ determines whether the active material is flow 
tumbling or flow aligning. The relaxation to equilibrium 
of Qij is determined by the molecular field Hij, which 
originates from the equilibrium free energy of an elastic 
anisotropic fluid, Hij = + (^*.^'/^)'^^5^7' where: 



T = L{^kQ^Jf/2 



(5) 



+ AQijQji/2 + BQijQjkQki/^ + C{QijQjk)'^ /A. 

L is the elastic constant, and A^B^C are material con- 
stants. Finally, the generalized advection term Sij and 
effective compressibility term P are defined as: 



P = P^p-L{^kQ^Jf/2, 



where Dij = {djUi ^ diUj)/2^ Qij = {djUi — diUj)/2 and 
Pq is a constant related to the hydrostatic pressure. For 
more details about the (non-active) model refer to j33|. 

We solve the dynamic equations for the active flow Ui 
and the orientation order Qij by using a hybrid Lattice 
Boltzmann (LB) algorithm. The time-evolution of Qij 
is obtained using an explicit finite difference scheme in 
time, whereas the equation for Ui is solved by a D3Q15 
lattice Boltzmann method [3^ . The simulations are per- 
formed in a circular simulation slice through a circular 
capillary tube with a cubic lattice; typically 4 x 200 x 
200 mesh points. Along the capillary axis (x), we as- 
sume periodic boundary conditions for both Qij and Ui^ 
whereas at the capillary walls we take no-slip boundary 
conditions for Ui (standard bounce-back) and fixed (in- 
plane or homeotropic) boundary conditions for Qij. All 
calculations are performed at small Reynolds numbers Re 
(typically Re ~ 10~^). The results are presented in units 
of the nematic correlation length (f, the typical length 
scale over which the nematic order can vary, which allows 
for application of the results to different active systems. 
Typically, the correlation length ~ a few times the length 
of the elementary constituents of the nematic and con- 
trols the characteristic size of the topological defects. For 
example, in bacterial suspensions j27|, ^ ~ 10 //m. Un- 
less otherwise stated, the following values for the material 
parameters are used: L = 40 pN, flow aligining regime 
^ = 1, F = 7.29/Pas, A = -0.2L/^^ B = -2.3L/^^ 
C = 1.91//^^, mesh resolution A = 1.5^, time step 
At = 2.2^V(^7). and capillary radius R = 150^. These 
parameters and simulation regime would be appropriate 
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FIG. 1: Active flow in a capillary with homogeneous boundary 
conditions, (a) Flow regimes upon increasing activity: basic 
state with no flow, first excited state with unidirectional flow, 
and second excited state with bidirectional flow. (b,c) Pri- 
mary flow Ux across the capillary diameter; secondary flow in 
the yz plane; and director profile of the active nematic (blue - 
out of plane, and red in plane) for (b) the unidirectional flow 
at C = 0.05, and (c) bidirectional flow at C = O-l- 



for dense bacterial suspensions [275, studied in microfiu- 
idic channels. 

Figure [1] shows the active flow profile in a capillary with 
surfaces imposing homogeneous alignment, parallel to the 
capillary axis x. Upon increasing the activity (Fig. [TJi), 
we observe a transition to spontaneous flow. Above the 
transition (first panel of Fig.[T]3), the leading flow is along 
the axis of the capillary, {x or — x), depending on the ran- 
dom noise in the inital condition for the director (0.1%). 
Interestingly, at higher activities ( (Fig. [it), the flow 
evolves into a bidirectional pattern, with the velocity in 
opposite halves of the capillary in opposite directions. 
The unidirectional flow can be interpreted as a first ex- 
cited active mode and the bidirectional flow as a second 
excited mode, selected by stronger active forcing. In is 
interesting to note that the stability of the two regimes 
partially overlap and they can be observed at the same 
activity (. This indicates the existence of energy barriers 
between active flow states. Indeed, analogous flow pat- 
terns and metastability of flow regimes was also observed 
in simulations of two-dimensional active nematics (sHj . 

However the full three dimensional simulations also al- 
low us to identify an inplane, secondary flow field within 
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FIG. 2: Active +1/2 defects are stabilised by cylindrical con- 
finement in a capillary. Active flow (in blue) and director 
(in red) for (a) homeotropic (perpendicular) and (b) circular 
director alignment at the capillary walls. The director and 
flow lie fully in the (y, z) plane, (c) Position of the defects A 
as a function of the activity ^ for homeotropic and circular 
surface, (d) Velocities Uy and Uz for homeotropic and circular 
surface at C = 0.001 along the y and z axes, respectively. 



the yz capillary cross-section (Fig. [TJd and c, second pan- 
els). The secondary flow consists of a circularly symmet- 
ric pattern of distinct vortices, 8 or 12 for the unidirec- 
tional and bidirectional flows, respectively. The vortices 
are pair- wise counter-rotating, thus exhibiting no net an- 
gular momentum. The magnitude of the secondary flow 
is 1% of that of the primary flow. The complex struc- 
ture of the flow is reflected in the active nematic direc- 
tor, shown in Fig. [TJd and c, third panels. The director 
exhibits a variable splay deformation in the radial direc- 
tion of the capillary and breaks the cylindrical symmetry 
in the bidirectional flow. Interestingly, the observed pri- 
mary and secondary active flows are separately analogous 
to the flow phases predicted in active polar films [36|. 

Topological defects in the orientational profiles of ac- 
tive materials could be used as sources of the microscopic 
activity, and a cylindrical capillary can naturally sup- 
port such active defect profiles. Two active defect con- 
figurations are shown in capillaries with surfaces impos- 
ing homeotropic (Fig. [2^) or circular (Fig. [JJd) director 
alignment. In both configurations, the two defect regions 
appear as a result of the orientational order being frus- 
trated by the confining surfaces and are actually defect 
lines along the capillary with +1/2 winding number. Sur- 
rounding the defects, a strictly in plane (%^z) active flow 
develops with zero flow component along the capillary 
axis. The magnitude of the flow is largest in the region 



surrounding the defects, which indicates that the defects 
act as effective local 'pumps' for the flow. It is the strong 
elastic distortion and variation in the degree of order in 
the defect cores, which produces locally high stress and 
consequently strong active flow. The small structures in 
the velocity field at the position of the defects seen in 
Fig. [2]i occur because local changes in the magnitude of 
the order parameter also cause elastic and active stresses. 
Flow along the direction of the defects (z) is larger by a 
factor of 2 than the flow in the perpendicular direction 
(?/). Note that the orientation of the defects clearly deter- 
mines also the direction of the local active flow (compare 
Fig. [2^ and Fig. [2)3). And depending on this direction, 
the defects are advected by the flow either towards the 
centre or the walls of the capillary (Fig.[2t). If defects ap- 
proach the surface, they produce larger distortions in the 
nematic and thus generate stronger flows. Quantitative 
flow profiles along the y and z capillary axes are shown 
for the homeotropic and circular boundary anchoring at 
C = 0.001 in Fig. [2]i. Finally, although the calculations 
are performed in a full three-dimensional capillary, the 
specific configuration with defects in Fig. [2] can be in- 
terpreted as quasi-two-dimensional active nematic in a 
circular confinement, since both flow and director have a 
strictly zero component along the capillary axis. 

At large activities, the active flow and the active de- 
fects escape in the third dimension (along the capillary 
axis), as indicated in Fig. [St and demonstrated in Fig.O 
The director flips out of the yz plane, and thus accommo- 
dates the distortion imposed by the surface by escaping 
in the third dimension rather than with the two +1/2 
defects. The active flow simultaneously develops a pri- 
mary x-component along the capillary axis (Fig.[3t). For 
C = 0.001, the ratio between the maximum secondary 
flow iiJTt^^ and maximum x-component of the flow u!^^^ 
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^ ju^^^ = 0.013. The nematic order is distorted 
and therefore there is no transition back to a zero flow 
state or an in plane director configuration upon reducing 
the activity back to zero. Thus the escaped flow regime 
can also be stable at the same lower activities (C 0) 
where defect states were observed (Fig.[2j), which is an in- 
dication of (meta) stable director fields -and consequently 
flow- states with energy barriers. 

The escaped active flow has, as for the homogenous 
capillary, a complex in-plane secondary flow. A lattice 
of vortices forms, with the number of vortices depending 
on the activity: 8 vortices for C, = 0.001 (Fig. [3li) and 
16 vortices for ^ = 0.01 (Fig. [3^). Figure [3]F depicts the 
primary flow along x direction and the secondary flow in 
the yz plane by showing flow streamlines. This reveals 
helical flow motifs, i.e. a lattice of vortices with a distinct 
pitch. The pitch of the active vortices is po ^ 11000^ for 
C = 0.001. The helical flow is reminiscent of cytoplasmic 
streaming in plant cells [37]. 

We next present results to determine the scaling of the 
active flow with the physical size of the capillary. Fig- 
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100 200 300 400 500 600 
^ Capillary diameter 2R / ^ 

FIG. 3: Escape of the active flow in third dimension, (a) 
Escaped active nematic director (red is in-plane, blue out-of- 
the-plane). (b) Velocity and director angle profiles across the 
capillary, (c) Primary and (d,e) secondary active flow profiles 
escaped along the capillary axis, (f) Stream-lines showing full 
active flow reveal an array of helical vortices. For visualisa- 
tion, the pitch is rescaled by a factor of 50 and data is drawn 
for 100 simulation boxes, which were periodically multiplied 
along the capillary axis, (g) Maximum flow in the centre of 
the capillary at C = 0.01 as a function of the capillary radius. 



ure[3^ shows the maximum flow velocity i^^^^, in the cen- 
tre of the capillary and along the capillary axis x, for vari- 
ous diameters of the capillary 2R at a given = 0.01. The 
dependence u^^^(K) is to good approximation linear, in- 
dicating a predominantly Poiseuille-like flow regime for 
the escaped active flow. The maximum velocity of a 
general isotropic fluid in a circular channel scales as 
'^^^^liso = fR^/^V^ where / is an effective body force 
driving the flow, e.g. a pressure drop per unit length of 
a capillary / = Ap/L. Indeed, the active contribution in 
the generalised Navier-Stokes equation (|2]) and dH can be 
extracted as an effective force from the active stress ten- 



sor /, = djU-f-- = -CdjQ^j = -{C/R)djQ,j a {(/R), 
where assuming scalalability of Qij the scale independent 
coordinates Xi = Xi/R have been introduced. By writ- 
ing u^^^ oc {( / R) R'^ / 4:7] (X i?, the linear scaling of u'^^^ 
is obtained, in agreement with the simulation results in 
Fig. [3^. Importantly, the escaped active nematic orienta- 
tional profile is well scalable in space and thus gives linear 
y^max ^ which is not necessarily true for other active 
nematic geometries or profiles. Note also that the slope 
of the u^^^{R) curve is proportional to C/?7, and could 
be used to measure this ratio in active nematic materi- 
als. Indeed using an alternative viscosity measurement 
the strength of the activity could also be determined, 
and the full relation between u^^^ and ("/ry (not only the 
proportionality) could be found by measuring at differ- 
ent C/^ ratios, e.g. by changing the temperature, yet 
retaining the escaped nematic profile. 

As experiments on active material in microchannels 
are increasingly feasible, there are several possible appli- 
cations for our results. Flow patterns similar to those 
seen here may emerge in dense suspensions of swimmers 
confined to microfluidic channels which have a preferred 
orientation relative to the surface [38| . Dense suspensions 
of active colloids provide another possible relevant exper- 
imental system [39|. Moreover similar flow patterns are 
seen in systems of collectively moving microtubules (ioj . 

In conclusion, we have demonstrated the richness of 
active nematic flow and orientational profiles in a cylin- 
drical capillary tube. The physical behaviour depends 
strongly on the surface boundary conditions imposed on 
the director field: we considered homogenous planar, cir- 
cular planar, and homeotropic (perpendicular) anchor- 
ing. The active flow emerged above a threshold in the 
activity for intrinsically homogeneous orientational pro- 
files, but at (" = for intrinsically inhomogenous director 
fields. The primary active flow developed along the cap- 
illary, but with intricate secondary flows, consisting of 
an array of velocity vortices. At high activities, bidirec- 
tional primary flow was found, indicating multiple ac- 
tivity modes which are observed to be (met a) stable rel- 
ative one to another. Topological defects imposed by 
the boundary conditions in this extensile active nematic 
study were shown to act as effective local pumps for 
the flow, with their topological structure determining the 
pumping direction. We demonstrated the escape of the 
active flow and active defects into the third dimension, 
showcasing the importance of the dimensionality in ac- 
tive materials, and providing a possible route to measur- 
ing the activity of active suspensions. 
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Curie Programme ACTOIDS. 
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